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We report on lattice calculations of the _B meson decay constant and the B-parameters using the NRQCD action for heavy quark. 
We performed a set of calculations at three lattice spacings in the quenched approximation in order to investigate the systematic 
uncertainty associated with the lattice discretization. Our results indicate that the systematic error is under good control for fg 
and Bg. We also present preliminary results from our unquenched simulation. 



1 Introduction 

The mass difference AM, in the neutral B q (q—s or d) 
meson system may be used to determine the CKM ele- 
ments \V t d\ or \Vtsl, and could eventually be a probe to 
explore new physics beyond the Standard Model. Ex- 
perimental effort to measure these quantities has already 
reached excellent accuracy for AM d B, and AM S is ex- 
pected to be observed in the near future. In order to 
determine \Vt q \ and then to search for the new physics 
effects using the mass difference, one has to know the 
matrix element of the four-quark operator 



(1) 



with B" and B» mesons as external states. Convention- 
ally, the matrix element is parameterized in terms of the 
lcptonic decay constant /r and the i?-parameter B B as 



(B q \O Lq ^ b )\B° q ) = y Bq B Bq (^ b )M Bq , 



(2) 



when the operator Oi, is defined at scale \i b . The mass 
difference is, then, given bya 



AM Bq = IC ■ \V t lV tq \ 2 [a s (fi b )} 
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a s (fJ, b 
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(3) 



with K, a known factor. The scale \x b is usually set to 
around m b . 

In this report we present recent results from a lattice 
QCD calculation of the B meson leptonic decay constant 
and the B-parameter obtained by the JLQCD collabo- 
rajjgjijgA part of this study has already been published 

It is also possible to estimate the width difference of 
B s mesons using the heavy quark expansion and lattice 
QCD. Our attempt is also reported at this conference!]. 



2 Lattice simulations 

To obtain reliable results from lattice calculations, sys- 
tematic uncertainty has to be under control. To disen- 
tangle possible systematic errors, we perform three sets 
of lattice simulations, which we call "NfO-A" , "NfO-B" 
and "Nf2" , respectively. The lattice parameters for these 
three sets are summarized in Table |l|. In these runs the 
physical volume of the lattice is kept at about (2 fm) 3 . 
Run "NfO-A" consists of three quenched lattices with dif- 
ferent lattice spacings, which may be used to study scal- 
ing behavior of the decay constant and the .B-parameter. 
Runs "NfO-B" and "Nf2" , both of which are still prelim- 
inary, have similar parameters except for the number of 
dynamical flavors. We use these runs to investigate the 
sea quark effect. 

In this work we use the lattice NRQCD action El in- 
cluding entire 1/m corrections for heavy quark. This 
approach has an advantage that the b quark can be di- 
rectly simulated and then extrapolation in heavy quark 
mass is not necessary. The light quark_is described by 
the 0(a)-improved Wilson quark actionB and, in partic- 
ular, in Runs "NfO-B" and "Nf2" the 0(a) -improvement 
is done nonperturbatively. 

Since the light quark much lighter than the strange 
quark can not be treated on the lattice with reasonable 
computer time, the simulations are carried out at four 
or five different light quark masses around and above 
the strange quark mass, and the results are extrapolated 
to the physical light quark mass. In particular, in the 
dynamical quark simulations "Nf2", independent gener- 
ations of gauge configurations are necessary to simulate 
different sea quark masses. 

3 Decay constant 

3.1 quenched results 

Several years ago, we performed a lattice calculation of 
the B meson decay constant using the NRQCD action! 3 ! 
The calculations in the quenched approximation were 
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Run 


NfO-A 


NfO-B 


Nf2 


N f 
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P 


6.1 


5.9 


5.7 


6.0 


5.2 


lattice size 


24 3 x64 


16 3 x48 


12 3 x32 


20 3 x48 


20 3 x48 


a [fm] 


0.086 


0.120 


0.176 


0.108 


0.095 


# of conf. 


518 


419 


420 


180 


3000 traj. 



Table 1: Simulation parameters. 




a [GeV 



Figure 1: Lattice spacing dependence of the decay constant in the 
quenched approximation. Error bars include statistical and sys- 
tematic errors except for the quenching. Open circles are from our 
previous work. Open triangle is a new result with nonperturba- 
tively 0(a)-improved light quarks (Run "NfO-B"). 



carried out at three lattice spacings to investigate the 
scaling property. The results are plotted in Figure [j] 
(open circles). 

At the present work, in Run "NfO-B" we performed 
another quenched calculation, in which the light quark 
action is nonperturbatively 0(a)-improved in contrast to 
the previous study, in which the improvement is done 
at the one-loop level. The result is also plotted in Fig- 
ure Q (open triangle). We observe from this figure that 
within the given uncertainty results are consistent with 
each other. It indicates that the systematic uncertainty 
depending on the lattice spacing is well controlled with 
the use of the 0(a)-improved actions. 

There are many other calculations of the B meson 
decay constant in the quenched approximation, for which 
recent summaries are presented in 113. It is remarkable 
that results from many groups with different approaches 
for heavy and light quarks are in reasonable agreement, 
including our results with the NRQCD. 



3.2 unquenching effects 

Last year, the JLQCD collaboration started a new 
project, namely QCD simulations with two flavors of dy- 
namical miarks. A preliminary report on this work is 
found inM. The calculation of the B meson decay con- 
stant is in progress on these simulations, which we call 
Run "Nf2" in this report. Unfortunately, at this stage 
we have not yet accumulated enough statistics to draw 
definite conclusion, so we present a status report of this 
work. All the results presented here are very preliminary. 
Because of the lack of statistics, the chiral limit of 



a s {M J 
a s (M Bd ) 



Pa) 



2//3o 



(4) 



still has large uncertainty as shown in Figure There 
we attempt various fit functions to obtain the physical 
light quark mass limit, and linear (solid) and quadratic 
(dotted) fits, are shown in Figure |[ Both of them have 
large x 2 /dof and the chiral limit is not very trusted. 
Perhaps, the reason for the bad statistical behavior of the 
unquenched calculation is that the simulation has long 
autocorrelation in the direction of fictitious time in the 
molecular dynamics evolution. More work is necessary 
to understand it and eventually to perform reliable chiral 
extrapolation. 

In spite of the poor data, especially toward the chiral 
limit, we attempt to see the effect of unquenching by com- 
paring the quenched and unquenched calculations of the 
decay constant. In Figure |^, the unquenched result (cir- 
cles) from Run "Nf2" are compared with the quenched 
one (triangles) from Run "NfO-B". For the unquenched 
data, we use the result of quadratic fit in the chiral ex- 
trapolation. It seems evident that the unquenched result 
is higher thanJie quenched one, as recently found by 
several groups More work (and statistics) is, how- 
ever, necessary to draw definite conclusion. Results with 
higher statistics will appear in the summer conferences 
in this year. 
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Figure 2: Chiral behavior of the heavy-light decay constant from 
Run "Nf2". The solid and dotted lines correspond to the linear 
and quadratic fits, respectively. 
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Figure 3: Comparison of the decay constant from quenched (tri- 
angles) and unquenched (circles) simulations. 
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Figure 5: BB d (jrib) as a function of 1/Mp at three lattice spac- 
ings (from Run "NfO-A") . The method 4 is used. The open and 
filled triangles in the static limit come from the lattice HQET ciaA- 
culations by Gimencz and Reyes with Nf=0 and 2, respectively 113. 



4-1 quenched results 

Quenched resultsJbr Bb on the Run "NfO-A" were al- 
ready reported in0. Figure |J shows the 1 /M dependence 
of BB d {rrib) obtained at ,9=5.9. Four symbols represent 
four different-methods to calculate the B-parameters as 
discussed in Q, and those have different systematic er- 
rors at order cCsAqcd/to or at order a 2 s . We may use 
the difference among the methods to estimate the size of 
systematic errors. 

In Figure |5| we plot all the results from Run "NfO- 
A" , where open circles, boxes and diamonds represent 
the results at /3=5.7, 5.9 and 6.1, respectively. 

We find that the results from different lattice spac- 
ings agree perfectly, suggesting that the systematic error 
associated with the lattice discretization is small. An- 
other result from Run "NfO-B" is also plotted in the same 
figure (open triangle-up), and we find a nice agreement 
again. 

The open triangle-down at 1/Mp—O is a result from 
the static approximation by Gimenez and Reyestd, which 
is consistent with our result extrapolated to the static 
limit. @ 
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Figure 4: Bs d {'mb) & s a function of 1/Mp at /3=5.9 in Run "NfO- 
A". Four symbols correspond to four different methods adapted. 



a We also note that the APEEJ and UKQCDL3 collaboration have 
recently calculated the £?-parameters using the relativistic lattice 
action for heavy quark, and the results are consistent with ours. 
Their approach has completely different systematic errors growing 
as a power of aM , heavy quark mass in the lattice unit. Therefore, 
an extrapolation is necessary from lighter heavy quark mass region 
to mi,. 
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4-2 unquenching effects 
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Using the same technique we developed in the quenched 
approximation, we are calculating the B-parameters in- 
cluding the effect of dynamical quarks in the Run "Nf2" . 
The problem of ill-behaved chiral extrapolation does not 
exist for Bb, since the light quark mass dependence of the 
-B-parameters is very small. Our preliminary result for 
BB d (mb) is shown in Figure^ (filled triangle-up) , where 
we do not find significant effect of the unquenching. 

Aj-urevious unquenched calculation by Gimenez and 
ReyesE3 is also plotted in the figure (filled triangle-down). 
There result seems significantly lower than others, but 
the calculation is done without the 0(a)-improvement 
and the result should have different systematic uncer- 
tainty. 

Our preliminary results for Bb are 



B B Mb) = 0.83(3)(8), 
B Bs /B Bd = 1.01(2)(2). 



(5) 
(6) 



5 Applications 



Let us briefly discuss the application of our results to the 
phenomenology. 

Substituting our results for Bs d (uLp) and the current 
world average for /s^by C. BernardE9 to the formula of 
mass difference Eq. (0), we obtain 



0.71 < 



v; b v td 



v; b v cd 



< 1.01. 



(7) 



Although the mass difference of B s mesons has not 
been discovered yet, precise measurement is expected in 
the near future at Tevatron Run II. We predict it as 

AM S = 18.2(4.7) ps" 1 . (8) 

The current experimental lower bound§ is 15.0 ps _1 . In 
the previous report we quoted higher value AM S — 
20.0(4.7) ps _1 . That is because the world average for 
fs, changed since that time. 

6 Conclusion 

It is essential for the determination of |Vtd| to calculate 
the B meson decay constant and the i?-parameters with 
controlled errors. The systematic errors associated with 
the lattice discretization are proved to be under good 
control in our calculation using the NRQCD action for 
heavy quark. Both for fg and Bb we show that the 
scaling violation is small in the quenched approximation. 

Moreover, we started a simulation including the ef- 
fect of dynamical quarks. At the moment, our statistics 
is not enough to draw definite conclusion, but we expect 
results in the near future. 
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